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ABSTRACT 

Galaxy redshift surveys are becoming increasingly important as a dark energy probe. We 
improve the forecasting of dark energy constraints from galaxy redshift surveys by using the 
"dewiggled" galaxy power spectrum, Pd w (k), in the Fisher matrix calculations. Since Pd w (k) 
is a good fit to real galaxy clustering data over most of the scale range of interest, our approach 
is more realistic compared to previous work in forecasting dark energy constraints from galaxy 
redshift surveys. We find that our new approach gives results in excellent agreement when 
compared to the results from the actual data analysis of the clustering of the Sloan Digital Sky 
Survey DR7 luminous red galaxies. We provide forecasts of the dark energy constraints from 
a plausible Stage IV galaxy redshift survey. 

Key words: cosmology: observations, distance scale, large-scale structure of universe 



1 INTRODUCTION 



One of the most important discoveries in m odern cosmology 
is the accelerated exp ansion of the Universe (Ries s et all 1 19981 : 
Perlmutter et al ]|l999l) . The power spectrum or 2-point correlation 
function measured from galaxy redshift surveys has provided one 
I ' of the primary probes of cosmic acceleration, both through the 
, broadband m easurement of the shape imprinted by matter-radiation 
■ equality (e.g. IPercival et al.112001 t lTegmark| [2004) and through the 
y—i ' baryon-a coustic oscillation (BA O) feature imprinted at recombina- 
, tion (e.g. lEisenstein et alj2005h . Galaxy clustering also allows us to 
. ' differentiate smooth dark energy and modified gravity as the cause 
for cosmic acceleration through the simultaneous measurements of 
!_h ■ the cosmic expansion history H (z), and the grow t h rate of cos- 
. "ff. 1 mic large scale s tructure, f g (z) dGuzzo et al.ll2008l : IWandl2008d : 
iBlake et alj2012h . 

The Fisher matrix approach has generally been used in the 
forecasts of future galaxy redshift surveys. In this paper, we im- 
prove the Fisher matrix approach by making it more realistic . This 
enables its use in cross-checking dark energy and gravity con- 
straints from current galaxy clustering data, as well as in making 
the forecasts for future galaxy redshift surveys more robust and re- 
liable. 

We present our method in Section [2] our results in Section [3] 
and summarize and conclude in Section|4] 
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2 METHOD 

The redshift-space galaxy power spectrum P{k,fj) is a rich 
source of cosmological information. It includes th e BAO feature 
jBlake & Gla zebrook 2003; Seo & Eisensteinll2003h . which has re- 
ceived a great deal of attention as a standard ruler that can be used 
in both the transverse direction (to measure distances) and the radial 
direction (to measure the Hubble rate). However, the full galaxy 
power spectrum at large scales is also sensitive to the underlying 
matter pow er spectrum, t o the growth of structure via redshift-space 
distortions dKaisenl 19871) . and to standard ruler effects. This addi- 
tional information requires some work to extract, since one must 
simultaneously measure the cosmology and the galaxy biasing pa- 
rameters. Nevertheless, the galaxy power spectrum provides the 
most powerful constraints on dark energy and gravity. In this paper, 
we focus on the analysis of the full set of 2-point galaxy statistics, 
and do not limit ourselves to only the BAO information. 



2.1 Formalism 

Our Fisher m at rix a pproach is de rived fr om t hat of 
ISeo & Eisensteml d2003l), an d based on IWand ( bOOfj. l2008j. 
l20ld) and lWangetal] feOlCh . In the limit where the length scale 
corresponding to the survey volume is much larger than the scale 
of any features in the observed galaxy power spectrum P g (k), we 
can assume that the likelihood funct ion for the band powers of a 
galax y redshift survey is Gaussian (Feldman, Kaiser. & Peacock] 
1 19941) . with a measurement error in In P(k) that is proportional to 
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\Vett (k)] 1//2 , with the effective volume of the survey defined as 



VeB (fc, /i) 



n(r)P g (k,n) 



i(r)P g (k,fi) + l 



nP g (k, n) 



nP g (k,fi) + 1 



Vsur 



(l) 



where the comoving number density n is assumed to only depend 
on the redshift (and constant in each redshift slice) for simplicity in 
the last part of the equation. 

In order to propagate the measurement error in In P g (k) into 
measurement er rors for the parameters Pi, we use the Fisher matrix 
dTegmarkll 19971) 



ainP 9 (k) 91nP 9 (k) 



dp t 



V cff (k) 



dk d 



2(2tt) 3 



(2) 



where pi are the parameters to be estimated from data, and the 
derivatives are evaluated at parameter values of the fiducial model. 
Note that the Fisher matrix Fij is the inverse of the covariance ma- 
trix of the parameters pi if the pi are Gaussian distributed. 

We adopt the standard notation that k can be decomposed into 
a line-of-sight component ku and the transverse or in-the-plane-of- 
the-sky component k± . The cosine of the angle between k and the 
line of sight vector is denoted by fi = feii/|k|. 



2.2 The model for the galaxy power spectrum 

At cosmological distances, the "true" galaxy power spectrum is not 
a direct observable, since one can measure a galaxy's position only 
in angular and redshift coordinates and not in its true 3D comoving 
coordinates. This is of course the basi s for extraction of the "stan- 
dard ruler" information, including the lAlcock & Paczvnskil ( 1 19791) 
effect. Therefore standard practice is to project the galaxies to their 
comoving positions assuming some reference cosmology (or fidu- 
cial cosmology), and then a power spectrum or correlation function 
estimator is applied. The observed galaxy power spectrum is then 
related to the true galaxy power spectrum via a coordinate transfor- 
mation: the wavenumber k rcf in the reference cosmology is related 
to the wavenumber in the true cosmology via 



^ k± and fcii cf = 



H lel (z) 



A 



(3) 



e (2) - ii H{z) IP 

Ba sed on ISeo & Eisensteinl d2003h and IChuang &~W ang 
d2012j) . our model for P 9 (k) can then be written as 



[D*f(z)] 2 H(z) + 



[D A (z)] 2 H^(z) 



xP dw (k) z e 



1 + k 2 /j, 2 al P /2 
*+-P sh ot, (4) 



where H(z) = a /a (with a denoting the cosmic scale factor) is 
the Hubble parameter, and Da(z) = r(z)/(l + z) is the angular 
diameter distance at z, with the comoving distance r(z) given by 



r(z) = c |S7fc | 1/2 si: 



,1/2 



dz 



H(z') 



(5) 



where sinn(x) = sin(:r), x, sinh(:r) for fife < 0, Qk = 0, and 
flfc > respectively. In addition to the geometrical distortion, this 
model includes the linear galaxy bias and redshift-space distortion 
(RSD), nonlinear smearing of the BAO feature, halo shot noise, 
small-scale peculiar velocities, and redshift errors. 



The bias between galaxy and matter distributions is denoted 
by b( z). The linear RSD parameter /3(z) = f g (z)/b(z) dKaised 
1 19871) . where L(z) is the linear growth rate; it is related to the lin- 
ear growth factor G(z) (normalized such that G(0) = 1) as follows 



/»(*) 



dlnG(z) 
din a 



(6) 



We have assume that the peculiar velocities of galaxies can be mod- 
eled with a probability distribution 



/(«) = 



<Tp\/2 



(7) 



where a v is the pairwise peculiar velocity dispersion. The Fourier 
transform of f(v) is + fc 2 ;x 2 cr 2 „/2 ], the small scale RSD fac- 
tor included in Eq. © dHamiltonll 19981) . Note that a r ,p is the dis- 
tance dispersion corresponding to the physical velocity dispersion 
Op, thus a p = H(z)[a(z)a r , P ], and 

(8) 



H(z)a(z) 



Note that we have adopted minimal small scale RSD modeling in 
this work (see Eq.Q), since we only consider quasilinear scales 
for a conservative app roach. The limitatio ns of Eq.((7} have been 
discussed in detail by IScoccimarrol (|2004). When smaller scales 
(< 20/i -1 Mpc) are included in the a nalysis, it will be critical to 
use an improved RSD model, see, e.g., IChuang & Wangl d2012bl) . 

An additional damping factor, e M CTr z , is inserted to ac- 
count for redshift uncertainties, with oy jZ = (dr/dz)a z . This is 
intended to incorporate the true redshift uncertainty resulting from 
fitting the centroid of an emission line (in an emission line survey), 
but this factor could also absorb other small errors in the redshift 
(e.g. due to the emission line velocity not being exactly equal to 
zero in the rest frame of the galaxy's host halo). 

Nonlinear smearing of the BAO feature occurs due to the 
small-scale (i.e. -C sbao ~ 150 Mpc) displacements during struc- 
ture formation. These displacements take sharp, coherent features 
in the correlation function at large scales (e.g. the BAO) and smear 
them out; in Fourier space, this corresponds to a damping of the 
oscillatory part of P(k), This effect is modeled by using the dewig- 
gled matter power spectrum at redshift z, given by 



Pdw(M) = G 2 {z)P k" s Ti w (k,z). 
Here T* w (k, z) is given by 

T d 2 w (k, z ) = r 2 (fc) e - ff " fe2/(2fc * ) 

+T 2 w (fc) [l - e -^ 2 /(^)" ( 



(9) 



(10) 



where T(k) is the linear matter transfer functio n, T nv/ (k) is the 
pure C DM (no baryons) transfer function given by Eisenstein & Hu 
( 19981 Eq. 29), and 



5M (k, z) = G 2 (z){l - M 2 + M 2 [l + / 9 ( 2 )] 2 } 



(11) 



describes the enhanced damping along the line of sight due to 
the enhanced power. The nonlinear dam ping factor, e~ gfjfc ^ 2k *\ 
with g u given by Eq.dl It. was derived bv lEisenstein. Seo. & White] 
(2007) using N-body simulations. Note that since density pertur- 
bations grow with cosmic time, the linear regime expands as we 
go to higher redshifts. Hence the function g M scales with the linear 
growth factor G(z) squared, which corresponds to the scale of the 
linear regime increasing with 1/G(z) at high redshifts. 

The scale fc« is rela ted to the percentage of nonlinearity from 
ISeo & Eisensteinl d2007l) . pnl, via 
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_1 = 8.355 /i _1 Mpc (er 8 /0.8) Pnl- 



(12) 



The true galaxy power spectrum should have pn l = 1. Recently 



"reco nstruction" algorith ms have been propos ed ( Eisenstein et al.l 
2007) and implemented JPadmanabhan et all |2012|) that reverse 
some of the flows and move galaxies back closer to their origi- 
nal (Lagrangian) positions. If such an algorithm is applied to data, 
the nonlinearity percentage can be reduced. BAO reconstruction is 
a rapidly developing field, but is in its early stages and high-z red- 
shift surveys may have to deal with survey geometries that are more 
complex and bias- weighted galaxy densities b 2 n that are smaller 
than that of e.g. BOSS. For the present work, we consider a range 
of values forpNL- The optimistic case of pnl = 0.5 corresponds to 
fc* ~ 0.24 /i/Mpc, whereas the most conservative case of pnl = 1 
(no reconstruction) corresponds to fc* ~ 0.12/i/Mpc, assuming 
as = 0.8. 

For an intuitive understanding of the dewiggled power spec- 
trum of Eq.((9}, we can rewrite its corresponding transfer function, 
Eo.dlOb. as follows 

- 9t ,fc 2 /(2fc2) 



Tdw(k, 



(*)+ [T 2 (fc)-T n 2 w (ft)] 

( fc)e -^ 2 /(^) 



(13) 



where we have defined T% AO (k) = T 2 {k) - T n 2 w (fc), the differ- 
ence between the linear matter transfer functions with and with- 
out baryons. Clearly, the exponential damping due to nonlinear ef- 
fects is only applied to the transfer function associated with BAO. 
Angul cTet alj J2008h have compared the spherically-averaged form 
of this model with measurements from simulated data, and found 
that it works extremely well on the linear and quasilinear scales; the 
assessment of its accuracy is presently limited by the shot noise of 
currently available numerical simulations. In future work, we will 
test this model without spherical averaging using numerical simu- 
lations to fully assess it. We do not expect this model to continue 
working well on the smallest scales, where the nonlinear damp- 
ing is coupled with RSD (e.g. [l ennings, Baugh, & Pascoli; [201ll ; 
iReid & Whitd201ltlChuang & Wandl2012bl) . 

To avoid the direct measurement of the unknown galaxy bias 
b(z), w e rewrite our model for the measured galaxy power spec- 
trum as dWandl2012h 



= P g (fc r x ef ,feJ e/ )/(/ l - 1 Mpc) 3 

[D A (z) Tct ] 2 H(z) 2 
- [<t b {z) + f g {z)a m {z) n \ 



[D A (z)] 2 H(z 
k 



Mpc 



where we have defined 

61/2 



Tdw(k, z) 



+ P S 1 



l + fcW, P /2 " shot ' 
(14) 



a g {z) = b{z)G{z)P Q L ^ and a m (z) = G(z) Pp 



1/2 



(15) 



The dimensionless power spectrum normalization constant Po is 
just Pq in Eq. l|9j in appropriate units: 



Po 



~ (Mpc/h) 3 (Mpc)™= I h n *' ( } 

The second part of Eq. d!6t is relevant if as is used to normalize 
the power spectrum. Note that 



> _ Tn s +2 

dk — -^-T 2 (k- /iMpc" 1 ) 



3ji(8fc) 



Xk 



(17) 



where k = k/[hMpc 1 ], and ji(kr) is spherical Bessel function. 
Note that Io = Io(u} m ,Ub, n 3 , h). Since ku and k± scale as H(z) 
and 1/ D A (z) respectively, Pg bs (k) in Eq. d 14t does not depend on 
h. 

Eg.ll 14t is the model we will use in this pap er. Its absorption 
of the bias factor is analogous to the approach of iSong & Percivall 
(2009), who proposed the use of f g (z)os(z) to probe growth 
of large scale structure. The difference is that Eg, (1141 uses 
f g (z)am(z) = f g (z)G(z)P 1 ^ 2 , which does not introduce an ex- 
plicit dependence on h (as in the case of using f g (z)as{z)). 



2.3 Parameters and assumptions 

In our method, the full set of parameters that describe the 
observed P 9 (k) are: {lnP(zi), InDA(zi), ln[f g (zi)a m (zi)}, 
lna g (zi), P s ' hot ; U) m , uj b , n a ,k,,a z /(l + z)}, where i indi- 



cates the i-th redshift slice, and cj„ 



Q, m h , and cut = 



Q b h 2 . We marginalize over {\na g (zi), P* hot } in each red- 
shift slice, as well as fc» and a z /(l + z), to obtain a Fisher 
matrix for {In H(z l ), In D A (zi), ln[f g (z l )a m (z l )]; u m , uj b , n s }. 
This full Fisher matrix, or a smaller set marginalized over various 
parameters, is projected into the standard set of cosmological pa- 
rameters {wo,w a , Qx, uj m , tot, n s , In A 3 }. Th ere are four d if- 
ferent ways of utilizing the information from P(k) (Wang 2012). 

It is important to note that when evaluating the derivatives of 
P g (k re ) with respect to the parameters described above (required 
to calculate the Fisher matrix), we should not extract information 
from the damping factors due to systematic uncertainties, in or- 
der to adhere to a conservative and robust approach. These damp- 
ing factors are only included to represent the loss of information at 
small scales due to nonlinear effects (and, if applicable, redshift un- 
certainties). We treat these damping factors as follows when deriva- 
tives are taken. 

The g M in the nonlinear damping factor, e _9f,fc '( 2fe *\ is 
fixed at fiducial model values when derivatives are taken, to avoid 
deriving cosmological information from the NL damping itself. 
Note that k is scaled as we vary Td w (k) for consistency, and we 
marginalize over fe* to allow for the significant uncertainty in the 
NL damping. 



The damping factor due to redshift uncertainty, e 



-k 11 CT^ 



computed with dr/dz from the fiducial model, to avoid deriving 
cosmological information from the damping itself. We marginalize 
over a z /(l + z) to allow for the uncertainty in our knowledge of 
redshift accuracy. 

The RSD factor due to small scale random motion of galaxies, 
1/[1 + k 2 n 2 a 2 iP /2], is fixed at fiducial model values when deriva- 
tives are taken. Note that as we vary H(z), this RSD factor remains 
unchanged, since k/i = fey oc H(z), while oy iP oc 1/H(z). The 
RSD factor is included here to represent the supression of power 
due to galaxy peculiar velocities, and not to provide an accurate 
modeling of RSD on all scales. 

Since our model fits real data well on these scales 
dChuang & Wandl2012al) . it represents a step forward in making 
Fisher matrix forecasting for galaxy redshift surveys more realis- 
tic. 



3 RESULTS 

We will present results on 

x h (z) = H(z) s/c and x d (z) = D A (z)/s, (18) 
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where s = r s (zd) is the sound horizon at the drag epoch, which is 
the characteristic scale of BAO. 

We assume the fiducial model adopted by the FoMSWG 
dAlbrecht et al. 2009), itself based on the 5-year Wilkinson Mi - 
crowave Anisotropy Probe (WMAP) results jDunklev et alj i2009): 
cj m ee tt m h 2 = 0.1326, ui b = tt b h 2 = 0.0227, h = 0.719, 
Q fe = 0, w = -1.0, n s = 0.963, and er 8 = 0.798. 

We will first present results for SDSS DR7 LRGs, in order to 
compare with the results from actual data analysis. The analysis of 
current GC data require the assu mption of cosmological p riors; we 
impose the same broad priors as Chuang & Wang] d2012al) . 

Next, we will present results for StagelV+BOSS spectro- 
scopic galaxy redshift surveys, and compare these with those 
from the previously w idely adopted approach derived from 
ISeo & Eisensteinl J2007b and developed in detail in IWand d2012h . 
No priors are used in deriving {xh{z), Xd(z), fg(z)o rn (z)/s a } 
constraints for future surveys, since these provide model- 
independent constraints on the cosmic expansion history and the 
growth rate of cosmic large scale structure. These allow the detec- 
tion of dark energy evolution, and the differentiation between an 
unknown energy component and modified gravity as the causes for 
the observed cosmic acceleration. 

In order t o derive dark energy figure of merit (FoM), as defined 
by the DETF dAlbrecht et alfcOOq) . we project our Fisher matrices 
into the standard set of dark energy and cosmological parameters: 
{u>o, w a , fix, fifc, Wmj ^b, n s , In A s }. To include Planck priorsQ 
we convert the Planck Fisher matrix for 44 parameters (including 
36 parameters that parametrize the dark energy equation of state in 
redshift bins) from the FoMSWG into a Planck Fisher matrix for 
this set of dark energy and cosmological parameters. 



3.1 Comparison with analysis of data 

To gauge the accuracy of our forecasting methodology compared to 
the full analysis of real data, we present our forecasts for the SDSS 
D R7 set of 87,000 LRGs in the redshift range 0.16-0.44 analyzed 
bv lChuang & Wand d2012al) in Tabled and compare them with the 
actual measurements performed as part of this work, using both the 
SDSS DR7 LRGs and the SDSS LR G mocks from LasDamafl 

The scale range analyzed by IChuang & Wand d2012al) is 
r =40-120 /i Mpc, w hich corresponds to the k = 2ir/r range of 
0.0524-0. 157 h Mpc" 1 . Ichuang & Wang] d2012ah used flat priors 
on ujb and n s that have widths of ±7ow map (with uwmap g iven 
by the W MAP seven year results from iKomatsu et al.l 4201 lb S. In 
addition, Chuang & Wang (2012a) imposed flat priors of 0.1 < 
j3 < 0.6, < a v < 500km/s, and 0.09 < k,(z = 
0.35)/[/iMpc" 1 ] < 0.13. We use Gaussian pr iors on the same 
parame ters with the same widths for the priors as Chuang & Wang 
d2012al) . and with means of a p = 250km/s, and &*/[/), Mpc -1 ] = 
0.11 G(z = 0.33) = 0.0939 (the width of fc»/[AMpc _1 ] is 
0.02 G(z = 0.33) = 0.01707). Note that our definition of k, is 
independent of redshift, thus it is divide d by the growth factor at t he 
effective redshift of the data set used by Chuang & Wang (2O12a)0 
In addition, we assume a redshift accuracy of a ln(l + z) — 5 x 
10" 4 , and a bias of b = 2.2 for the SDSS LRGs. These additional 



For a genera l and robust method for including Planck priors, see 
iMukheriee eta3 fe008l). 

2 URL: http : //l ss .ph y . vanderbilt . edu/lasdamas/ 

3 IChuang & Wand <2012iJ) scaled their results from z. ss = 0.33 to 
z e f f = 0.35 in order to compare with previous results by other groups. 



assumptions are not needed for the analysis of r eal or mock data, 
since t he overall amplitude is marginalized over (Chuang & Wang 
l2012ah . 

Since the Fisher matrix results depend on the fiducial model 
assumed, we give our Fisher matrix forecasts for two differ- 
ent fiducial models in Table Q] the FoMSWG fiducial model 
dAlbrecht etaUT 2009) with with u m = Q m h 2 = 0.1326, cu b = 
Q b h 2 = 0.0227, h = 0.719, n k = 0, w = -1.0, n s = 

0.963, and as 0.7 98, and the Euclid Red Book fiducial model 

dLaureiis et alj|20l"lh with with uj m = Q. m h 2 = 0.1225, cu b = 
Q. b h 2 = 0.021805, h = 0.7, n k = 0, w = -0.95, n s = 0.963, 
and as = 0.8. These fiducial models lead to a\nxh(z) and 
a In Xd{z) that differ by 3.1% and 6.8% respectively. 

The model used here , Eq.(ll4b . differs somewhat from that 
used by IChuang & Wang] d2012j). Our new model, E q.QD, 
uses aisotropic dewiggling whereas |chuang & Wand d2012a[) used 
isotropic dewiggling, which neglects the additional damping along 
the line of sight due to the enhanced Lagrangian displacement in 
redshift spaceQ Using our Fisher matrix method, we find that as- 
suming isotropic dewiggling leads to an under-estmate of cr\ nx . r z ) 
and er lna;£i ( z ) of ~ 23 — 24% and ~ 12 — 13% respectively. 

In order to make an accurate comparison, we have repeated the 
analysi s of the SDSS DR7 LRG sample used by 1 Chuang & Wangj 
d2012al) using Eq . d 1 4t as part of this work. We find that the data 
and mocks give o-\ nXh ( z ) and o~i nXd ( z ) that differ by 13% and 30% 
respectively, with the mock results agreeing with our Fisher matrix 
forecasts at a level of 10% or better, given the dependence of the 
Fisher matrix forecasts on the assumed fiducial model. 

Our Fisher matrix forecasts for the measurement uncertainty 
on f g (z)a m (z) /s 4 are in excellent agreement with the results from 
the MCMC analysis of the LasDamas SDSS LRG mocks, while 
the results from the MCMC analysis of SDSS DR7 LRG data give 
significantly smaller uncertainty on fg(z)a m (z)/s 4 . This is likely 
due to the apparent excess clustering of SDSS DR7 LRGs along 
the line of sight (this is apparent from compari ng the mock and 
data panels of Fig.l in Chuang & Wang 2012a), which is likely 
also responsible for the smaller than expected measurement uncer- 
tainty of Xh(z) and Xd(z) from this data sample. This excess power 
along the line of sight explains the widely noted excess power 
on large scales f or the spherically-average d galaxy cor r elatio n 
function (see, e.g ., Cabre & Gaztanagai (2009): K azin et alj f20101; 
IChuang, Wang. & Hemanth j d2012h ). Since the BOSS CMASS 
galaxies do not h ave the excess clustering on the same large scales 
dReid et al]|2012h . the excess large scale clustering of SDSS DR7 
LRGs must be due to sample variance or unknown systematic ef- 
fects. 

Taking all the factors discussed above into consideration, our 
Fisher matrix forecasts are in excellent agreement with the results 
from actual data analysis. It is reassuring that our Fisher matrix 
method gives very similar results compared to actual data analy- 
sis, making it a reliable tool for parameter forecasting for future 
surveys. 

3.2 Forecasts for a Stage IV galaxy redshift survey 

We perform forecasts for a Stage IV galaxy redshift survey cov- 
ering an area of 15,000 deg 2 , using slitless grism spectroscopy to 
detect the Ha emission line. A wavelength range of 1.1-2.0 /im, 
corresponding to 0.7 < z < 2.0, was assumed. The depth of the 

4 This is equivalent to setting g M — > G 2 (z) in Eq. lilt . 
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Model 


Method 


Xh{z) 







fg(z)cr m (z)/s 4 


FoMSWG 
EuclidRB 


Fisher matrix 
Fisher matrix 


731% 
7.09% 


4.99% 
4.67% 


21.98% 
22.61% 


20.69% 
21.05% 


None 


MCMC analysis of data 


5.80% 


3.74% 


14.89% 


14.01% 


None 


MCMC analysis of mocks 


6.64% 


5.37% 


23.72% 


22.61 % 



Table 1. Our Fisher matrix estimate of the percentage precision of measurement of 11,(2) = H(z)s/c, Xd(z) = Da(z)/s, /?, and f a (z )cr m (z) / at an 
effectiv e redshift of z = 0.35 from SDSS DR7 LRGs, compared to the actual measurements using the anisotropic correlation function per Chuang & Wang 
(2012a). Eq.fi4l is used in all cases. 



survey was computed using instrument parameters (throughput, ex- 
posure time, e t c.) sim ilar to those provided for the Euclid mission in 
lLaureiis et al.l d201ll) : it is thus representative of a next-generation 
space-based galaxy survey, though it may not correspond precisely 
to the final Euclid numbers. 

All our forecast results are shown for Stage IV plus BOSS. 
The BOSS survey is assumed to cover 10,000 (deg) 2 , a redshift 
range of 0.1 < z < 0.7, with a fixed galaxy number density of 
n — 3 x 10~ 4 /i 3 Mpc~ s , and a fixed linear bias of b = 1.7. 

We discuss the galaxy yields from a Stage IV galaxy redshift 
survey in detail in Sec l3.2.ll For clustering analysis, we also require 
the galaxy bias. We use the galax y bias function for emission line 
galaxies given bv lOrsi et alTd2010h . which increases with redshift 
reaching b — 1.7 at z = 2. Again, we note th at this is likely to be 
conservative: the recent bias determination of iGeach et al. I l2012h 
for Ha emitters is b — 2.4^2 at z — 2.23. We assume a redshift 
accuracy of a z /(l + z) = 0.001, and a peculiar velocity dispersion 
of <t p = 290km/s. 

We consider two different cutoffs in scale: k max = 0.2 fa/Mpc 
and kmax ~ 0.2 h/Mpc, in order to include the quasilinear regime 
only in our forecasts. The choice of k max = 0.2 /i/Mpc is con- 
servative, and represents the lower bound of the scale range in 
which our model works well in analyzing real data. The choice 
of kmax =0.3 ft/Mpc is more optimistic, but represents a feasible 
goal for the lower bound of the scale range in which future studies 
will enable robust and accurate modeling. 



3.2. 1 Galaxy yields for a Stage IV redshift survey 

Galaxy yiel ds were computed u sing the exposure time calculator 
described in iHirata et al.l ( T20121) . Two exposures on each field in 
each grism bandpass were assumed. The zodiacal background was 
set to that at 45° ecliptic latitude and 90° away from the Sun at the 
mean of the annual cycle, and we include a foreground dust column 
of E(B — V) = 0.05 magnitudes; these values vary over any real- 
istic survey but are representative. Standard read noise assumptions 
for the 2kx2k Teledyne HgCdTe detectors were used (32 channel 
readout, 1 frame per 1.3 s, 20 electrons rms per correlated double 
sample, with a noise floor of 5 electron rms for many reads). The 
galaxy survey was assumed to be 70 per cent complete down to 
the flux limit for a 7<r significance matched-filter detection^ The 
extinction-corrected Ha flux limit varies with redshift and galaxy 



size, but is in the range of (2.2 — 3.6) x 10 " erg s 'cm 
source half-light radius of 0.3 arcsec. 

As a point of comparison, we ran the lHirata et al. 1 J2012I) code 
on the Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) 
G141 grism (1.1 < A < 1.7 11m), using parameters from the In- 
strument Handbook ( lDresselll20T II). We find that for an exposure 
time of 2700 s, the 5a sensitivity should be (4.6 - 8.2) x 10~ 17 
erg s -1 cm -2 , with the lower (better) numbers at the red end of the 
bandpass. This is in good agreement (~ 20 per cent) with the me- 
dian sensiti vity actually achi eved by WFC3 observations - see e.g. 
Figure 5 of lAtek et al] d2010h . 

The line flux sensitivity and completeness are only part of 
determining the number of redshifts obtained by a survey - one 
also needs a luminosity function. In the past decade of space- 
based redshift survey mission planning, the Ha luminosity func- 
tion (HaLF) has been a matter of vigorous debate: direct measure- 
ments have suffered from small-number statistics, while indirect 
methods (based on scaling from rest-frame ultraviolet or [O II] 
luminosities) have had difficult-to-quantify systematic errors. In- 
deed, the estimates used fo r space mission planning (e.g. Yan et al] 
Il999l: iHopkins et al]|200fj|; iReddv et al]|2008l ; Ijouveletal] 12009} 
IGeach et al . 2010; Sobral et al. 2012) have spanned a factor of ~ 3 
in number density, even accounting for the different cosmologies 
assumed. Fortunately, empirical measures of the HaLF across the 
relevant range of redshifts with large-number statistics (dozens of 
objects in the relevant flux range, in multiple fields) are now avail- 
able. 

We us e the HaLF of Sobral etlill d2012l) . with conversions 
described in lHirataetal.ld2012h §3E to ensure consistency with the 
exposure time calculator inputs. This is based o n blind narrow-ban d 
surveys, and updates the previous estimate bv lGeach et al.ld2010h . 
The new HaLF is lower than the previous estimate; in approximate 
decreasing order of importance, the main differences are: 

1. Consistent treatment of internal (host galaxy) extinction correc- 
tions, which are applied to some HaLF results and must be undone 
to predict redshift survey yields. 

2. Improved statistics and addition of data at new redshifts. 

3. Redshift-averaging effects in some of the grism luminosity func- 
tions (this does not occur in narrowband surveys). 

4. Aperture corrections. 

5. Conversion to the WMAP-5/FoMSWG cosmology. 



5 Note that some forecasts in the literature use other definitions of detection 
significance, based on other extraction apertures. The differences are often 
tens of percents and occasionally as large as a factor of 2. The matched-filter 
method gives the highest reported significance. 



The narrowband surveys do not cleanly separate the Ha 6563 A 
line fr om the [N II] doublet at 6548,6583 A, and the Sobr al et all 
d2012h HaLF removes the estimated [N II] contribution. Of course, 
in a grism survey the two lines will be a partial blend, thus we may 
be underestimating the final detection significance of the galaxies. 
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For this reason, we expect that our analysis is somewhat conserva- 
tive. 

Table[2]gives our resultant galaxy yields as a function of red- 
shift. 



3.2.2 Dark energy figure of merit results 

Table [3] shows the dark energy figure-of-merit (FoM) dWand 
20083), 



FoM{wa,w a ) 



-y/det Cov{w , w a ) 



(19) 



for the four different approaches to utilizing the i nformation from 
the measured anisotropic galaxy power spectrum dWandl2012l) : 
a){x h (z),x d (z)}fromP(k); 

(2) {x h (z),x d (z), f g (z)a m (z)/s a } from P(k); 

(3) P(k), marginalized over f g (z)a m (z); 

(4) P(k)+/ S (z); P(k) including f g (z)a m (z). 

It is clear from Table [3] that for a given cutoff kmax, the 
FoM for (wo,w a ) increases as we increase the dewiggling scale 
ht, (i.e., decrease the nonlinear effects). For a fixed level of non- 
linearity (i.e., fixed fc*), the FoM for (wo,w a ) increases as we 
increase the cutoff k ma x. The scaling of f g (z)a m (z) with s de- 
pends on the level of nonlinearity assumed: for 50% nonlinear- 
ity (fe» = 0.24/i/Mpc), a ~ 4, while for 100% nonlinearity 
(fc* = 0.12 /i/Mpc), a ~ 5. This is not surprising, since the scal- 
ing of f g (z)am{z) with s 4 (i.e., a = 4 ) originates from the lin- 
ear matter power spectrum dWandl2012h . When nonlinear effects 
are fully included (and not assumed to be reduced due to density 
field reconstruction), the appropriate model for P(k) (i.e., Eg. 1141 ) 
deviates significantly from the linear power spectrum, leading to 
modification of the scaling of f g (z)a m (z) with s. 

The choice of a does not affect the FoM(wo,w a ) 
from the {x n , xj, j g o m j s a } case, as the correlations between 
f g (z)(j m (z) I s a and {xh(z),Xd(z)} depend on a. Choosing the 
a that minimizes the uncertainties in f g (z)a m ,(z)/s a does maxi- 
mize the FoM(«7o, w a ) when Planck priors are included. 



3.3 Comparison with previous work 

Previously, the forecasts of dark energy constraints from full P(k) 
assumed that 



P° ld {k r l f ,k\ ef ) 



[Da(z)] 2 H(z)™f 
xb 2 (l + p M 2 ) 2 P Un (k|z) 



xe 2 



(20) 



where the linear matter power spectrum Pu„(k\z) — 
G 2 (z)P k" s T 2 (k) (with T{k) denoting the linear matter 
transfer function), and 



2 

(7 



(I) H?)' 

Alternatively, we can write 



(21) 



pold/uref U re f\ 
r g (n ± ) 



P 9 oid (fcl e/ ,fc; e/ )/(^ _1 Mpc) 3 



[Da(z 



i re/1 



■H(z) 



[o-g(z) + f g {z)a m {z) ^1 



l^Mpc" 

+ Pshot- 



T\k)e 



(22) 



Table |4] lists the FoM for (wo, w a ) for the same four cases as 
listed in Table [3] Each line in Tableland its corresponding line in 
Table [3] assume the same level of nonlinearity and the same cutoff 
k m ax ■ The only difference between the two tables is the model as- 
sumed for P(k): Eq. d 14b (from Eq. |4)) is assumed for Table [3] 
while Eq. \22\ (from Eq. 1201 ) is assumed for Table|4] 

Note that the two assumed models of P(k) give similar FoM 
for all the cases that marginalize over the growth information, and 
for the cases that include growth information but assume only a 
nonlinearity level of 50% (pnl = 0.5 or k* = 0.24 h/Mpc). 
When we assume a nonlinearity level of 100% (pnl = 1 or 
fc* = 0.12/i/Mpc), our new model (Eq. 11141 ) gives significantly 
larger FoM for the cases that include the growth information. This 
is because the old model in Eq. d22t simply damps the linear mat- 
ter power spectrum exponentially, while the new model in Eq. d 1 4b 
only damps the BAO oscillations, and retain smaller scale infor- 
mation via the "no- wiggle" matter power spectrum P n „(k|z) = 
G 2 (z)P k ns T 2 w (k) (wifhT^Cfc) denoti ng the zero baryon trans- 



Eisenstein & Hul dl998l) ). This results in sianif- 



fer function from 
icantly smaller uncertainties in ln/3(z) (and In f g (z)a m (z) / s a ) 
when our new model is used, which in turn leads to significantly 
larger FoM(^o, w a ) when growth information is included. 

While the scaling of f g (z)a m (z) with s depends on the level 
of nonlinearity assumed in our new model (see discussion in the 
previous sub section), f g (z)a m (z) scales with s 4 in the old model 
dWandl2012l) . We find that our FoM results are not sensitive to the 
exact choice of a; we have chosen a — 4 for all the FoM tabulated 
in Tables[3]and|4]when growth information is included. 



4 CONCLUSION 

We have shown that the forecasting of dark energy constraints from 
galaxy redshift surveys can be improved in fidelity by using the 
"dewiggled" galaxy power spectrum, Pd w (k), in the Fisher matrix 
calculations. Since Pj w (k) is a good fit to real galaxy clustering 
data over most of the scale range of interest, our approach is more 
realistic compared to previous work in forecasting dark energy con- 
straints from galaxy redshift surveys. 

We tested our methodology by comparing our Fisher matrix 
forecasts with results from actual data analysis, and found excel- 
lent agreement (see Table QJ. Our Fisher matrix method gives very 
similar results compared to actual data analysis, making it a reliable 
tool for parameter forecasting for future surveys. 

Using our new approach, we studied a Stage IV galaxy red- 
shift survey, in combination with BOSS, without and with Planck 
priors. We find that in this new approach, increasing nonlinear ef- 
fects from 50% (best case) to 100% (most conservative) has a sig- 
nificantly reduced impact on the dark energy figure of merit com- 
pared to previous work. This indicates that the erasure of informa- 
tion by nonlinear smearing is having only a modest effect on our 
ability to constrain cosmology using the "full P(k)" method in our 
new realistic approach. 
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z A EE50 dVKdz-dA) F Hm @0.30" n dNIdz-dA 

(fim) (arcsec) (Mpc 3 deg -2 ) (Wm -2 ) (Mpc~ 3 ) (deg~ 2 ) 



0.700 


1 


.1160 


0.2297 


5.56334E+06 


3.56250E- 


19 


1.70876E- 


-04 


9.50640E+02 


0.750 


1 


.1489 


0.2315 


6.06057E+06 


3.33879E- 


19 


1.89962E- 


-04 


1.15128E+03 


0.800 


1 


.1817 


0.2335 


6.54313E+06 


3.15425E- 


19 


2.11061E- 


-04 


1.38100E+03 


0.850 


1 


.2145 


0.2354 


7.00894E+06 


3.01278E- 


19 


2.23002E- 


-04 


1.56301E+03 


0.900 


1 


.2474 


0.2373 


7.45644E+06 


2.89615E- 


19 


2.01417E- 


-04 


1.50185E+03 


0.950 


1 


.2802 


0.2393 


7.8845 1E+06 


2.79942E- 


19 


1.81900E- 


-04 


1.43419E+03 


1.000 


1 


.3130 


0.2413 


8.29240E+06 


2.72837E- 


19 


1.62666E- 


-04 


1.34889E+03 


1.050 


1 


.3458 


0.2433 


8.67969E+06 


2.67418E- 


19 


1. 4493 IE- 


-04 


1.25796E+03 


1.100 


1 


.3787 


0.2454 


9.04621E+06 


2.63267E- 


19 


1.28950E- 


-04 


1.16651E+03 


1.150 


1 


.41 15 


0.2474 


9 39200F+06 


2 60703F— 


19 


1.13930E- 


-04 


1 07003F+03 


1.200 


1 


.4443 


0.2494 


9.71731E+06 


2.64718E- 


19 


9.39210E- 


-05 


9.12659E+02 


1.250 


1 


.4771 


0.2696 


1 00225F+07 

I .\J\J ±- j— ivy / 


2 80452F— 


19 


6 69603F- 


-05 


6 71 109F+02 


1.300 


1 


.5100 


0.2714 


1.03081E+07 


2.68660E- 


19 


6.78461E- 


-05 


6.99362E+02 


1.350 


1 


.5428 


0.2731 


1.05746E+07 


2.59173E- 


19 


6.78039E- 


-05 


7.16996E+02 


1.400 


1 


.5756 


0.2749 


1.08226E+07 


2.51215E- 


19 


6.73077E- 


-05 


7.28446E+02 


1.450 


1 


.6084 


0.2768 


1.10529E+07 


2.44917E- 


19 


6.61038E- 


-05 


7.30641E+02 


1.500 


1 


.6413 


0.2786 


1.12662E+07 


2.39935E- 


19 


6.34722E- 


-05 


7.15093E+02 


1.550 


1 


.6741 


0.2804 


1.14632E+07 


2.35608E- 


19 


6.03252E- 


-05 


6.91521E+02 


1.600 


1 


.7069 


0.2823 


1.16446E+07 


2.31880E- 


19 


5.73065E- 


-05 


6.673 13E+02 


1.650 


1 


.7397 


0.2841 


1.18112E+07 


2.29092E- 


19 


5.40911E- 


-05 


6.3888 1E+02 


1.700 


1 


.7726 


0.2860 


1.19637E+07 


2.26675E- 


19 


5.11258E- 


-05 


6.1165 1E+02 


1.750 


1 


.8054 


0.2879 


1.21027E+07 


2.24586E- 


19 


4.83888E- 


-05 


5.85637E+02 


1.800 


1 


.8382 


0.2898 


1.22291E+07 


2.22811E- 


19 


4.58519E- 


-05 


5.60729E+02 


1.850 


1 


.8710 


0.2917 


1.23435E+07 


2.21390E- 


19 


4.34464E- 


-05 


5.3628 1E+02 


1.900 


1 


.9039 


0.2936 


1.24465E+07 


2.20333E- 


19 


4.11589E- 


-05 


5.12286E+02 


1.950 


1 


.9367 


0.2955 


1.25388E+07 


2.20230E- 


19 


3.86029E- 


-05 


4.84035E+02 


2.000 


1 


.9695 


0.2975 


1.26210E+07 


2.20856E- 


19 


3.59559E- 


-05 


4.53799E+02 



Table 2. Galaxy yields for a 2-exposure Stage IV galaxy redshift survey as discussed in Sec. 13.2.1] Columns indicate: the redshift; observer-frame wavelength 
of Ha; the half-light radius (encircled energy 50%, EE50) of the point-spread function; the cosmological volume element dV/ (dz ■ dA); the limiting flux for 
a 0.3 arcsec half-light radius galaxy; the number n of observed sources per unit comoving volume; and the number of sources dN/dz-dA per unit redshift per 
unit solid angle. 







FoM 


FoM Gfl (FoM, d7) 


FoM 


FoMq r (FoM, d7) 


(TiMpc- 1 ) 


(hMpc" 1 ) 






{^(k)} 




{P(k)+f g } 


0.2 


0.12 


6.56 


30.78 (22.80, 0.0514) 


14.81 


40.48 


(24.47, 0.0476) 


0.2 


0.24 


10.05 


45.30 (31.07, 0.0470) 


23.37 


54.59 


(35.05, 0.0456) 


0.3 


0.12 


9.83 


44.11 (33.98, 0.0437) 


19.94 


65.32 


(35.19, 0.0386) 


0.3 


0.24 


12.73 


62.51 (42.41, 0.0394) 


29.81 


79.57 


(46.43, 0.0374) 


Stage IV+BOSS+Planck 


0.2 


0.12 


58.30 


139.24 (80.49, 0.0392) 


61.30 


171.90 


(106.98, 0.0341) 


0.2 


0.24 


92.64 


193.61 (119.58, 0.0376) 


96.22 


238.63 


(152.10, 0.0329) 


0.3 


0.12 


85.24 


209.98 (110.10, 0.0344) 


89.29 


240.11 


(136.75, 0.0311) 


0.3 


0.24 


119.88 


273.20 (152.31, 0.0322) 


123.55 


315.23 


(184.21, 0.0295) 



Table 3. Ou r Fisher matri x forecasts for Stage IV+BOSS galaxy redshift surveys using our new galaxy power spectrum model, Eq. dl4t . for the four cases 
discussed in Wang (2012). FoMg^j denoted the FoM assuming general relativity. The parameter 7 is defined by fg(z) = [Sl m (a)] 7 . 
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